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Abstract 

We explore the phase diagram and the critical behavior of QCD thermodynamic quantities in 
the context of the so-called Polyakov-Nambu-Jona-Lasinio model. We show that this improved 
field theoretical model is a successful candidate for studying the equation of state and the critical 
behavior around the critical end point. We argue that a convenient choice of the model parameters 
is crucial to get the correct description of isentropic trajectories. The effects of the regularization 
procedure in several thermodynamic quantities is also analyzed. The results are compared with 
simple thermodynamic expectations and lattice data. 

PACS numbers: 11.30.Rd, 11.55.Fv, 14.40.Aq 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) exhibits at zero temperature and chemical potential 
two remarkable features which play an essential role in our understanding of strong interac- 
tion phenomena: the fundamental degrees of freedom are colorless bound states of hadrons 
(quark confinement), and chiral symmetry is dynamically broken. As it is well known these 
features characterize the nonperturbative nature of the QCD vacuum. It is expected that, at 
large energy densities, the so-called QCD phase transition occurs: the interaction becomes 
weaker and weaker due to asymptotic freedom with the formation of a new state of 

matter, the quark gluon plasma (QGP), and chiral symmetry gets restored. 

The study of the QCD phase diagram in the (T, /x)-plane and the search for signatures 
of the QGP have attracted an intensive investigation over the last decades. The outputs of 
this research are expected to play an important role in our understanding of the evolution 
of the early universe and of the physics of heavy-ion collisions at the Brookhaven National 
Laboratory, and at LHC (CERN), in the future. 

Various results from QCD inspired models indicate (see e.g. Refs. {2I, Q]) that at low 
temperatures the transition may be first order for large values of the chemical potential; on 
the contrary a crossover is found for small chemical potential and large temperature. This 
suggests that the first order transition line may end when the temperature increases, the 
phase diagram thus exhibiting a (second order) critical endpoint (CEP) 5|-|8( that can be 
detected [si, by a new generation of experiments with relativistic nuclei, as the CBM 
experiment (FAIR) at GSI. Fodor and Katz ll|] claim the values T CEP = 162 MeV and 
fj, = 360 MeV for such a critical point, although its precise location is still a matter 



of debate 



12J. In the chiral limit it is found, in accordance with universality arguments, a 



tricritical point (TCP) in the phase diagram, separating the second order transition line from 
the first order one. So, the exploration of the critical region of the phase diagram of strongly 
interacting matter gains increasing attention, both experimentally and theoretically. 

As an approach complementary to first-principle lattice simulations, one can consider 
several effective models. One of them is the Nambu-Jona-Lasinio (NJL) model, that is 
undoubtedly a useful tool for understanding chiral symmetry breaking but does not possess 
a confinement mechanism. As a reliable model that can treat both the chiral and the 
deconfinement phase transitions, we can consider the Polyakov loop extended NJL (PNJL) 



2 



model 



13 



16] . In the PNJL model the deconfinement phase transition is described by the 



Polyakov loop. This improved field theoretical model is fundamental for interpreting the 
lattice QCD results and extrapolating into regions not yet accessible to lattice simulations. 

A nontrivial question in NJL type models is the choice of the parameter set and of the 
regularization procedure. In fact, one should keep in mind that this type of models are 
used not only to describe physical observables in the vacuum but also to explore the physics 
at finite temperature and chemical potential. As it is well known, the order of the phase 
transition on the axis of the (T, /i)-plane is sensitive to the values of the parameters, most 
notably to the value of the ultraviolet cutoff needed to regularize the integrals. In the pure 
NJL model a large cutoff leads to a second order transition, a small cutoff to a first order one 
~y\ . An interesting feature to be noticed is that the requirement of global accordance with 
physical spectrum is obtained with values of the cutoff for which the transition is first order 
on the T = axis and a smooth crossover on the \x = axis of the phase diagram. However, 
it has also been shown that different parameter sets, although providing a reasonable fit 
of hadronic vacuum observables and predicting a first order phase transition, will lead to 
different physical scenarios at finite T and \x 18|, ll9[. For instance, the absolute stability of 
the vacuum state at T = is not always insured. 

The consequences of the choice of the parameter set for the scenarios in the (T, /x)-plane 
have not been discussed in the framework of the PNJL model, where the most popular 
parameter set does not allow for the absolute stability of the vacuum at T = 0. In the 
present work, our main goal is to analyze this problem and we will assume that the most 
reliable parametrization of both NJL and PNJL models positively predicts the existence of 
the CEP in the phase diagram, together with the formation of stable quark droplets in the 
vacuum state at T = 0. 

Finally, the physical relevance of our numerical solutions is insured by requiring the 
agreement with general thermodynamic requirements. In particular, we will verify that the 
correct description of isentropic lines is closely related with the parameter choice in the pure 
NJL sector. 

Concerning the regularization of some integrals, since, as it has been noticed by several 
authors, the three dimensional cutoff is only necessary at zero temperature, the dropping 
of this cutoff at finite T is carefully analyzed in this work: this procedure allows for the 
presence of high momentum quark states, leading to interesting physical consequences, as 
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it has been shown in [20|, where the advantages and drawbacks of this regularization have 
been discussed. We will enlarge the use of this procedure to the PNJL model and discuss 
its influence on the behavior of several relevant observables. 

Let us notice that the choice of a regularization procedure is a part of the effective 
modeling of QCD thermodynamic. Indeed the presence of high momentum quarks (no 
cutoff on the temperature dependent part of the loop integrals) is required to ensure that 
the entropy scales as T 3 at high temperature. Hence we found that a comprehensive study 
of the differences between the two regularization procedures (with and without cutoff on the 
quark momentum states at finite temperature) is necessary to have a better understanding 
of the PNJL model and the role of high momentum quarks around the phase transition. 
This is one of the main purposes of this paper. 

This paper is organized as follows: In Sec. II we present the model and formalism starting 
with the deduction of the self-consistent equations. We also extract the equations of state and 
the response functions, and show the pertinence and physical relevance of a convenient choice 
of the parametrization and regularization procedures of the model. Section III is devoted 
to the study of the phase transition at zero temperature, showing the important role of the 
choice of parameters for the formation of quark droplets in mechanical equilibrium with the 
vacuum at zero pressure. In Sec. IV we study thermodynamical quantities that, compared 
with lattice results, point out the physically relevant regularization procedure at T ^ 0. The 
enlargement to \i ^ allows for the study of the phase diagram in the (T, /i)-plane (Sec. 
V). In Sec. VI we discuss the important role of the choice of the model parameters for the 
correct description of isentropic trajectories. In Sec. VII we proceed to study the size of 
the critical region around the CEP and its consequences for the susceptibilities and critical 
exponents. Finally, some concluding remarks are presented in Sec. VIII. 

II. MODEL AND FORMALISM 

A. Model Lagrangian and gap equations 

The generalized Lagrangian of the quark model for Nf = 2 light quarks and N c = 3 
color degrees of freedom, where the quarks are coupled to a (spatially constant) temporal 
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a 


Ol 


«2 


^3 


^3 


bi 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



TABLE I. Parameter set used for the Polyakov loop potential (2) and (3). 



background gauge field (represented in term of Polyakov loops) , is given by 




Cpnjl = q(i^ l D^ - m)q + -g[(qq) 2 + (qi^rq) 2 ] 

-U(9[A],9[A]',T), (1) 

where the quark fields q = (u, d) are defined in Dirac and color spaces, and m = 
diag(m u , ma) is the current quark mass matrix. The pure NJL sector contains three parame- 
ters: the coupling constant g, the cutoff A, and the current quark mass m = m u = m^, to be 
determined (see Sec. II C) by fitting the experimental values of several physical quantities. 

The quarks are coupled to the gauge sector via the covariant derivative = <9 M — iA^. 
The strong coupling constant gstrong has been absorbed in the definition of A^{x) = 
9 strong A^(x)^f, where A% is the SU C (3) gauge field and A a are the Gell-Mann matrices. 
Besides in the Polyakov gauge and at finite temperature A^ = 5qA° = —i5%A. 

The Polyakov loop $ (the order parameter of Z 3 symmetric/broken phase transition in 
pure gauge) is the trace of the Polyakov line defined by: $ = -^-((Vexpi J dr A 4 (x,r) ))„, 
where ((■ ■ with /? = 1/T is the thermal expectation value in the grand canonical ensem- 
ble. 

The pure gauge sector is described by an effective potential U (&[A], 3? [A]; T) that takes 
the form 



W($,$;T) 



.MQ$$-|($» + $3) + ^($$) 2 



(2) 



where 



= a + ai ( ^ 



T 



T 



a 2 — + a 3 — 



Tr 



T 



(3) 



The coefficients To, a, and foj of the Polyakov loop effective potential are chosen (see Table 
[Hand |l5j) to reproduce, at the mean-field level, the results obtained in pure gauge lattice 
calculations. 
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From the Lagrangian ([T]) in the mean-field approximation it is straightforward (see 
Ref. 22|) to obtain the constituent quark mass M that is given by 



M = m — 2g (qq) , (4) 

where the quark condensate (qq) has to be determined in a self-consistent way. So, taking 
already into account Eq. (j3J), the PNJL grand potential density in the mean-field approach 
is given by 
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23|: 



/j3 
(Ztt) 

d 3 p 



+ 2N f T 



(2tt 



d 3 
(27T) 

, [\nN+(E p ) + In N$ (E, 



pj\ ' 



(5) 



where E p is the quasiparticle energy for the quarks, E p = \/ p 2 + M 2 , and by defining 



E p n, the upper sign applying for fermions and the lower sign for antifermions, the 



functions N£ and N$ are: 



n£(e p ) 

N*(E P ) 



1 + 3 + $ e -^ P + ) e -Pzi + e - 3 ^ + 
1 + 3 ($ + $>e- pE A e^^ + e~ 3/3 ^ 



(6) 
(7) 



It was shown in Ref. 



that all calculations in the NJL model can be generalized to the 
PNJL one by introducing the modified Fermi-Dirac distribution functions for particles and 
antiparticles: 



f+(E p ) = M+iV+, M+ = U + 2i>e-^ e"^ + + e" 3 ^ 



(8) 
(9) 



To obtain the mean-field equations we must search for the minima of the thermodynamical 
potential density ([5]) with respect to (qq), $ and $. In fact, by minimizing Q with respect 
to (qq), we obtain the equations for the quark condensate: 

d 3 p M 



(qq) 



-6 



■[0(A 2 



V 



ft(E p )-U(E p )\. 



(10) 



Let us stress that in the latter we used a cutoff A only for the T = part of the integral 
(denoted case I in the following); one could put a global cutoff J Q A (denoted case II in the 
above). 



Furthermore, minimization of Q with respect to $ and $ provides, respectively, the two 



additional mean-field equations 
for case I and J Q A for case II): 



23] (the integral that appears has to be understood as J Q C 



rp4 

~ T 

- 12T 



[-6 2 (T)$-6 3 $ 2 + 64*$ 2 ] 
d 3 p 



(27TY 



-PE- N - 



(11) 



~ T 

- 12T 



[-6 2 (T)$ - & 3 $ 2 + h$<S> 2 ] 
d 3 p 



(2vr) £ 



€ -Pe; N + + e -WE» N - 



(12) 



The solutions of the three coupled equations (jlj), ffTTT) and (|12[) allow us to obtain the 
behavior of M, and the Polyakov loop expectation values as a function of T and \i. 



B. Equations of state and response functions 

From the thermodynamical potential Q(T,fi) one can derive equations of state which 
allow us to compare some of our results with observables that have become accessible in 
lattice QCD at nonzero chemical potential. The relevant observables are the (scaled) quark 
number density, defined as 



p q (T,fj) i fdn\ 



T 



T 3 T 3 V dfi J , 

and the (scaled) "pressure difference" given by 

Ap(T,n) _ p(T,tj,)-p(T,0) 



(14) 



As usual, the pressure, p, is defined such as its value is zero in the vacuum state [18] and, 
since the system is uniform, one has 

P( T,,) = -^, (15) 

where V is the volume of the system. 

Our work also includes the study of the isentropic trajectories due to their relevance for 
the study of the thermodynamics of matter created in relativistic heavy-ion collisions. The 
equation of state for the entropy density, s, is given by 

s(T,p) = (%) , (16) 



dT 



and the energy density, e, comes from the following fundamental relation of thermodynamics 



e(T,/i) = Ts(T,p) + pp q (T,p) - p(T,p). (17) 

The energy density, as well as the pressure, is defined such that its value is zero in the 
vacuum state fis| . 

The baryon number susceptibility Xg an d the specific heat C are the response of the 
baryon number density p q (T, p) and the entropy density s(T, p) to an infinitesimal variation 
of the quark chemical potential p and temperature, given respectively by: 

dp q \ _ J ^ T fds\ 



*'=U?J T ' a,ld c= v\m) t - (18) 

These second order derivatives of the pressure are relevant quantities to discuss phase 
transitions, mainly the second order ones. 



C. Model parameters and regularization schemes 



As already stated, the pure NJL sector involves three parameters: the coupling constant 
g, the current quark mass m and the cutoff A. These parameters are determined in the 
vacuum by fitting the experimental values of several physical quantities. We notice that the 
parameters g and A are correlated with each other: if we increase g in order to provide a 
more significative attraction between quarks, we must also increase the cutoff A in order to 
insure a good agreement with experimental results. In addition, the value of the cutoff itself 
does have some impact as far as the medium effects in the limit T = are concerned. 

We remember that different parameterizations may give rise to different physical scenarios 
at T = and p ^ 18|, even if they give reasonable fits to hadronic vacuum observables 
and predict a first order phase transition. Here, we will use two different sets of parameters 
whose values are presented in Table [TTJ These two sets are the most widely used in NJL 



type models: set A is taken from 



6| anc 



in the context of the PNJL model 



15 



set B from [24J, the last one being commonly used 



231 ] . The main feature is a lower (larger) value of 



the cutoff for set A (B), for which we verify that A/M < 1.8 (> 1.8). We notice that the 
transition between the regime of stable to the regime of metastable quark matter occurs at 
the value A/M m 1.8 25j. So, we will prove that only the set of parameters A insures the 
stability conditions and, consequently, the compatibility with thermodynamic expectations. 
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A 
[GeV] 


9 

[GeV~ 2 ] 


m 
[MeV] 


|<VwM 1/3 

[MeV] 


fn 

[MeV] 


run 
[MeV] 


M 
[MeV] 


Set A 


0.590 


7.0 


6.0 


241.5 


92.6 


140.2 


400 


Set B 


0.651 


5.04 


5.5 


251 


92.3 


139.3 


335 



TABLE II. Set of parameters (A, g, m) used in the NJL sector of the PNJL model and the physical 
quantities chosen to fix the parameters. The constituent quark mass obtained is also included. 



We notice that set A also agrees with an empirical relation derived in 26] which states that 
stable quark matter is only possible in NJL model if M > 4f n . 

On the other hand, the regularization procedure, as soon as the temperature effects 
are considered, has relevant consequences on the behavior of physical observables, namely 
on the chiral condensates and the meson masses j^. In PNJL model, the two types of 



15 



23 



27 



regularization may be found (15J, |23|, |28H30| as well as different sets of parameters 
29j. So, in order to compare the differences between the use of different sets of parameters 
and of regularization schemes, in the physical scenarios of the PNJL model, we will consider 
sets A and B of parameters and two different regularization procedures at T ^ 0: 

Case I. — The cutoff is used only in the integrals that are divergent (A — > oo in the 



convergent ones; see Eqs. OlUJ [TT] fP2|) for example) at finite temperature, a 



procedure that 



15 



17 



2s 



3l|. 



23 



29]. 
20]. 



allows us to take into account the effects of high momentum quarks 

Case II. — The regularization consists in the use of the cutoff A in all integrals 
Advantages and drawbacks of these regularization procedures have been discussed in 
Here, our main goal is to show nontrivial consequences of the regularization scheme used in 
case I. We remind that the main drawback of this regularization is that at high temperature 
there is a too fast decrease of the quark masses that become lower than their current values. 
This leads to a non physical behavior of the quark condensates that, after vanishing at the 
point where constituent and current quark masses coincide, changes sign and acquires a 
nonzero value again. Therefore, if we want to keep calculating observables in this region, it 
seems sensible to impose the condition that the quark masses take their current values and 
the quark condensates remain zero. This is the approach used here. 
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III. PHASE TRANSITION AT ZERO TEMPERATURE 



Our study refers mainly to the finite temperature case. However, the particular case 
of zero temperature is very important due to the possibility of having, simultaneously, a 
vanishing density and a finite chemical potential. This feature depends on the choice of the 
parameters and is necessary in order to insure the satisfaction of general thermodynamic 
requirements. 

Here we analyze the stability of the quark matter at T = 0. For this special case the 



PNJL model reduces to the NJL one. We wi 
of the system along the same lines of 



18 



19 



1 now present a discussion about the stability 



32J, which will allow us to choose the set of 



parameters corresponding to the most convenient physical situation. 

In the limit T — > the normal Fermi-Dirac distribution function reduces to the step 
function 

9(ji - E p ) = 6(p F -p)6(p - M), (19) 
where the Fermi momentum is given by 



Vf = vV - M 2 6(p - M) = (tt 2 Pq ) 1/3 6(n - M). (20) 

The important point of our argumentation about the choice of the model parameters 
comes from the comparison between the point (0, p c ) of the phase diagram, where p c is the 
position of the first order line at zero temperature, and the point (0, M vac ), where M vac = M 
is the mass of the u, <i-quark in the vacuum. Two special cases are observed 



(i) For set A, the first order phase transition occurs at /i c such that p c < M vac , and 
consequently (see Eq. [20J) the phase transition connects the vacuum state (p q = 0) 
directly with the phase of partially restored chiral symmetry (p q = p c ). 

(ii) For set B, p c > M vac , so the phase transition connects a p q ^ phase of massive 
quarks with the phase of partially restored chiral symmetry (p q = p c ). 

So, although we can choose several sets of parameters which fit physical observables in the 
vacuum, we notice, however, that the value of the cutoff itself does have some impact on the 
characteristic of the first order phase transition. Comparing the two sets of parameters we 
verify that for larger values of the cutoff, as in set B, a more strong attraction is necessary 
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both to reproduce the physical values in the vacuum and to insure a first order phase 
transition. As we will argue in the sequel, the more reliable case is provided by set A. 
In case (i) the energy per particle reaches, at p q = p c , an absolute minimum e < e(0) = 



M vac . This is compatible with the existence of stable quark matter, indicating the possibi 
for finite droplets to be in mechanical equilibrium with the vacuum at zero pressure 



31 



18 



ity 



19|, 



33] . This is due to the fact that the pressure has three zeros, respectively at p = 
0, 0.52p , 4.3po (po = 0.17 fm~ 3 ), that correspond to extrema of the energy per particle. 
The third zero of the pressure, located at p c = 4.3p , corresponds to an absolute minimum 
of the energy. The critical point of the phase transition in these conditions satisfies to 



Pc < M, 



vac 



18 



34[ | . This can be seen by comparing p c = 383 MeV with the quark masses 



M vac = M = 400 MeV. Above p q = p c , we have again a uniform gas phase. For densities 
< p q < p c the equilibrium configuration is a mixed phase, where the equality of all 
intensive variables (T, P and p) defines the condition for the phase equilibrium. So, the 



Gibbs criterion is satisfied and the phase transition is a first order one 18 



On the contrary, the minimum at p q = p c in case (ii) corresponds to metastable quark 
matter, since e > e(0) = M vac . The pressure still has three zeros, respectively at p q = 
0, 1.51po; 2.64p , that correspond to extrema of the energy per particle. The main difference 
now is that the absolute minimum of the energy per particle is at p q = 0. This means that, 
in spite of being in the presence of a first order phase transition, there are no droplets in 
mechanical equilibrium with the vacuum at zero pressure. In addition, from the physical 
point of view, this scenario is unrealistic because it predicts the existence of a low-density 
phase of homogeneously distributed constituent quarks [3]. Other implications of this 
scenario on the reliability of isentropic trajectories will be discussed later. 



IV. THERMODYNAMIC QUANTITIES 

A significative information on the phase structure of QCD at high temperature is obtained 
from lattice calculations in the limit of vanishing quark chemical potential. The transition 
to the phase characteristic of this regime is related with chiral and deconfinement transitions 
which are the main features of our model calculation. 

Following the argumentation presented in [jjj, we use the reduced temperature T c by 
rescaling the parameter Tq from 270 to 190 MeV (we do this rescaling only for the remainder 
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in 




200 300 400 

T (MeV) 



500 600 



FIG. 1. Scaled pressure (p), energy per particle (e), and entropy (s) as a function of the temperature 
at zero chemical potential for both sets of parameters A and B, and both regularization procedures. 

of this section). In this case we loose the perfect coincidence of the chiral and deconfinement 
transitions: they are shifted relative to each other by less than 35 (30) MeV for set A (B). 



15], we define T c as the average of the two transition temperatures: we have 



As in Ref. 

T c = 190 (184) MeV for set A (B) within the range expected from lattice calculations [35]. 



For comparison purposes with lattice findings, we start by considering our numerical 
results at vanishing quark chemical potential by checking the usefulness of the present reg- 
ularization procedure, case I. To this purpose, we plot the scaled pressure, the energy and 
the entropy as functions of the temperature in Fig. [TJ 

The transition to the high temperature phase is a rapid crossover rather than a phase 
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transition and, consequently, the pressure, the entropy and the energy densities are contin- 
uous functions of the temperature. For case I we observe a similar behavior in the three 
curves: a sharp increase in the vicinity of the transition temperature and then a tendency 
to saturate at the corresponding ideal gas limit. Asymptotically, the QCD pressure for Nf 
massless quarks and — 1) massless gluons is given, for \x = 0, by 

^? = (Nl-l)^ + N c N f 7 -^, (21) 
T 4 45 180 V ; 

where the first term denotes the gluonic contribution and the second term the fermionic one. 
Our results follow the expected tendency and go to the free gas values, a feature that 



was also found with this type of re gula rization in the context of the PNJL model 



13 



For what concerns the NJL model 



15|- 



171 ] let us notice that if indeed a tendency to saturate 



is found, the asymptotic value is at about half the ideal gas limit. Hence the inclusion of 
the Polyakov loop effective potential $) (it can be seen as an effective pressure term 
mimicking the gluonic degrees of freedom of QCD) is required to get the correct limit. 

The inclusion of the Polyakov loop together with the regularization procedure imple- 
mented in case I, is essential to obtain the required increase of extensive thermodynamic 
quantities, insuring the convergence to the Stefan-Boltzmann (SB) limit of QCD 



Some 



comments are in order concerning the role of the regularization procedure for T > T c . In 
this temperature range, due to the presence of high momentum quarks, the physical situ- 
ation is dominated by the significative decrease of the constituent quark masses by the qq 
interactions. This allows for an ideal gas behavior of almost massless quarks with the correct 
number of degrees of freedom. 

The advantage of our phenomenological model is the possibility to provide equations of 
state at nonzero chemical potential too. So, we can also test its ability to reproduce recent 
progress in lattice QCD with small non vanishing chemical potential. 

In order to do this, we plot in Fig. 2 the scaled pressure and the quark number density 



for jj, = 0.6 T c (it was done in 15J for set B). The parameters are the same as in Fig. 1 and 



only case I for the cutoff procedure is considered. The agreement with the lattice data is 
fairly good, showing that the general pattern of these quantities and the behavior at large T 
is reproduced. The more significative deviations are observed in the scaled baryon number 
density at intermediate temperatures as is evident from Fig. 2 (right side). However, there 
are some advantage with respect to the employed parametrization A, where an improvement 
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0.8 



0.6 



0.4' 



0.2 



0.0- 



Case I: A — > c° 



H = 0.6T C : 

Set A 

— SetB 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 



T/T 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.( 

T/T c 



FIG. 2. Comparison of the scaled pressure difference (left panel) and the scaled quark number 
density (right panel), as a function of temperature at finite chemical potential for sets A and B, 



with the lattice data taken from Ref. 



46|. 



of the results is observed as shown in the left panel of Fig. 2. 

In conclusion, by introducing finite chemical potential we are able to compare the obtained 
results with lattice data and check the validity of both sets of parameters. From Fig. [2] we 
conclude that, for case I, both sets of parameters are in good agreement with lattice results, 
in particular, to Ap/T A and p q /T 3 at /i = 0.6 T c . 



V. PHASE STRUCTURE 



The phase diagram for both sets of parameters (see Fig. [3]) is determined by the behavior 
of the orderlike parameters (qq), $ and $ together with the grand canonical potential as 
a function of temperature and chemical potential. To draw the phase diagram we will use 
T = 270 MeV as given by pure gauge lattice calculations, a choice that ensures the very 
important physical outcomes of lattice calculations that chiral and deconfinement transitions 
coincide in the PNJL model.. 

We start our analysis in the limit T = 0, where the first order phase transition occurs 
at the same chemical potential for both cases I and II: fi = 383 (344) MeV for set A (B). 
This is due to the fact that at T = all the integrals are regularized with the cutoff A. 
As a matter of fact, in this limit the Fermi functions in the gap equations become a step 
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Parameter 
set 


Regularization 
procedure 


rpCEP 

[MeV] 


^CEP 

[MeV] 


T TCP [MeV] 
(chiral limit) 


fi TCP [MeV] 
(chiral limit) 


Set A 


Case I: A — > oo 


172.48 


286.35 


206.50 


192.87 




Case II: A = const 


169.11 


321.32 


207.66 


270.80 


Set B 


Case I: A — > oo 


87.99 


328.84 


162.39 


253.06 




Case //: A = const 


87.71 


329.51 


163.06 


268.05 



TABLE III. Location of the CEP and the TCP at the (T, /x)-plane for both sets of parameters and 
regularization procedures. 

function of the form 9 (A — pp) where pp is the hadronic matter Fermi momentum. So, the 
integration occurs between pp and A for both cases. 

From Fig. 3 we also see that, at p — 0, the crossover takes place for set A at T = 235 (272) 
MeV for case I (II), while for set B the crossover takes place at T = 229 (256) MeV for case 
I (II). 

At nonzero chemical potential, as the temperature increases, it is well know that the first 
order transition persists up to the CEP. At the CEP the chiral transition becomes of second 
order. For temperatures above the CEP a smooth crossover takes place. 

These general characteristics are qualitatively similar for both cases, in the two sets of 
parameters, as it was expected. The relevant point is the distance between the CEP's, in 
cases I and II, is bigger for set A than for set B. This is due to fact that the CEP's for set A are 
at a higher temperature than for set B. As the temperature increases, the high momentum 
quarks, that are taken into account in case I (A — > oo), are more and more relevant leading 
to a visible splitting of the lines of first order phase transition and, consequently, to other 
location of the CEP in the phase diagram. For set B this splitting is smaller, once lower 
critical temperatures are observed. In the chiral limit (m = and m n = 0), the transition 
is of second order at p — and, as p increases, the line of second order phase transition will 
end in a first order line at the TCP. The location of the tricritical points are also included 
in Table IIHI Nevertheless, in this case, the TCP is located at higher temperature than the 
CEP (see Fig. [31 right panel), and we already see a bigger shift between them. 
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FIG. 3. Phase diagrams for cases I and II in the chiral limit (boxes) or not (circles) using the 
parameters set A (left panel) and B (right panel). The solid part of the curves denotes a first 
order transition, the dashed part to the second order transition and the dotted line the crossover 
transition. 

VI. NERNST PRINCIPLE AND ISENTROPIC TRAJECTORIES 



The isentropic lines play a very important role in the understanding of thermodynamic 
properties of matter created in relativistic heavy ion collisions. Most of the studies on this 
topic have been done on lattice calculations for two flavor QCD at finite /x but there 
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Some model calculations predict 



are also studies using different type of models |34|, 
that in a region around the CEP the properties of matter are only slowly modified as the 
collision energy is changed, as a consequence of the attractor character of the CEP jo]. 

Our numerical results for the isentropic lines in the (T, /i) -plane are shown in Fig. [4j 
where we have used set A of parameters and both regularization procedures. 

We start the discussion by analyzing the behavior of the isentropic lines in the limit 
T — > 0. We point out that, as already referred by other authors js^j], in this limit: 

(i) s — > 0, according to the third law of thermodynamics; and 

(ii) for s/p q = const, we have to insure that also p q — > 0. 

However, the satisfaction of the condition (ii) is only provided when \i < M vac . In spite of 
the general use of set B in the literature of the PNJL model, only set A satisfies this ansatz. 



1(3 
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H (MeV) n (MeV) 

FIG. 4. Isentropic trajectories in the (T, /i)-plane for case I (left panel) and case II (right panel) 
using the parameter set A. The following values of the entropy per baryon number have been 
considered: s/p q = 1, 2, 3, 4, 5, 6, 8, 10, 15 (anticlockwise direction). 

We remember (Sec. 2.3) that with set A we are, at T = 0, in the presence of droplets (states 
in mechanical equilibrium with the vacuum state at P = 0). 

Consequently, even without reheating in the mixed phase as verified in the "zigzag" shape 



of 37H40|. all isentropic trajectories directly terminate at the first order transition line at 
T = 0. So, for set A it is verified that s — > and p q — > in the limit T — > 0, as it should be. 

In conclusion, our convenient choice of the model parameters allows a first order phase 
transition that is stronger than in other treatments of the NJL (PNJL) model. This choice is 

nn 

crucial to obtain important results: the criterium of stability of the quark droplets [18|, [19| 
is fulfilled, and, in addition, simple thermodynamic expectations in the limit T — y are 
verified. 

At T 0, in the first order 
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ine, the behavior we find is somewhat different from those 
41| where a phenomena of focusing of trajectories towards the 
CEP is observed. For case I (see Fig. HI left panel) we see that the isentropic lines with 
s/ p q = 1, ...,4 come from the region of symmetry partially restored and attain directly the 
phase transition, going along with the phase transition as T decreases until it reaches T = 0. 
The same behavior is found for case II when s/p q = 1,2 (see Fig. HJ right panel). For case 
II, we also observe, in a small range of s/p q around 3, a tendency to convergence of these 
isentropic lines towards the CEP. These lines come from the region of symmetry partially 
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restored in the direction of the crossover line. For smaller values of s/p q , the isentropic lines 
turn about the CEP and then attain the first order transition line. For larger values of s/p q 
the isentropic trajectories approach the CEP by the region where the chiral symmetry is 
still broken, and also attain the first order transition line after bending toward the critical 



point. As already pointed out in [34[, this is a natural result in these type of quark models 
with no change in the number of degrees of freedom of the system in the two phases. As the 
temperature decreases a first order phase transition occurs, the latent heat increases and 
the formation of the mixed phase is thermodynamically favored. 

In the crossover region, for both cases, the behavior of the isentropic lines is qualitatively 



similar to the one obtained in lattice calculations |37J or in some models 38|, |42|, |43| . On 
the other hand, the isentropic trajectories in the phase diagram indicate that the slope of 
the trajectories goes to large values for large T. We can also conclude that, in the PNJL 
model, the entropy and the baryon number density are very sensitive to the regularization 



procedure used 17J, |20[, and this effect is also relevant for the present situation. 



VII. SUSCEPTIBILITIES AND CRITICAL EXPONENTS 

The grand canonical potential (or the pressure) contains the relevant information on ther- 
modynamic bulk properties of a medium. Susceptibilities, being second order derivatives of 
the pressure in both chemical potential and temperature directions, are related to fluctua- 
tions that are supposed to represent signatures of phase transitions of strongly interacting 
matter. In particular, the quark number susceptibilities play a role in the calculation of 
event-by-event fluctuations of conserved quantities such as net baryon number. Across the 
quark hadron phase transition they are expected to become large, what can be interpreted 
as an indication for a critical behavior. We also remember the important role of the second 
derivative of the pressure for second order points like the CEP. 



In previous works 
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32], we have studied the CEP within the restrictions imposed by 
the regularization in case II. It is important to investigate if the type of regularization plays 
a significant role in the critical properties of physical observables, such as the baryon number 
susceptibility and the specific heat, and respective critical exponents, in the vicinity of the 
CEP. The relevance of these physical observables is due to the size of the critical region 
around the CEP which can be found by calculating the baryon number susceptibility, the 
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FIG. 5. Phase diagram in the PNJL (case I and case II) for the parameter set A. The size of the 
critical region around the CEP is plotted for XqlXq ee = 2,3,5. 

specific heat and their critical behaviors. The size of this critical region is important for 
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future searches of the CEP in heavy ion-collisions 

In our calculations, we will use only the set A of parameters due to the advantages of this 
set as explained in the previous sections. In Fig. [5] we plot the phase diagram in a region 
around the CEP for both cases. 

The way to estimate the critical region around the CEP is to calculate the dimensionless 
ratio Xq/Xq Tee y where x{ ree is obtained taking the chiral limit m — 0. Figure [5] shows a 
contour plot for three fixed ratios (Xq/x{ ree — 2.0, 3.0, 5.0) in the phase diagram around 
the CEP. We notice an elongation, in the direction parallel to the first order transition line, 
of the region where Xq is enhanced, indicating that the critical region is heavily stretched in 
that direction. 

The elongation of the critical region in the (T, /i)-plane, along the line of the phase 
transition (or the crossover), is larger in case I (see for example Xq/x{ ree = 2-0 in Fig. |5]). It 
means that the divergence of the correlation length at the CEP affects the phase diagram 
quite far from the CEP, particularly in case I, and a careful analysis including effects beyond 



the mean-field needs to be done 



44]. 



We remember that one of the main effects of the Polyakov loop is to shorten the temper- 



ature range where the crossover occurs 23[ . On the other hand, this behavior is boosted by 
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Quantity critical exponents/path 


oasc i 


Pooo TT 


Universality 




e 1 from left (=>■) 


0.66 ±0.01 


0.66 ±0.01 


2/3 


Aq 


e / from right (^=) 


0.69 ± 0.01 


0.69 ± 0.02 


2/3 




7 g / from left (— Y) 


0.53 ± 0.02 


0.51 ± 0.01 


1 /o 

1/2 






0.64 ±0.01 


0.62 ± 0.02 


2/3 




a. 1 from below (ff) 










a x = 0.58 ±0.01 


a x = 0.52 ±0.01 




c 


a' / from above (JJ-) 


0.68 ± 0.01 


0.68 ± 0.01 


2/3 




a 1 from below (t) 


0.50 ±0.01 


0.47 ±0.02 


1/2 



TABLE IV. The arrow =>(t) indicates the path in the £t(T)- direction to the CEP/TCP for 
fi < fi CEP (T < T TCP ). 

using A — > oo: at fi = the crossover occurs between about 180 MeV and 270 MeV, for case 
I, and between about 210 MeV and 325 MeV, for case II. The combination of both effects 
results in higher baryonic susceptibilities even far from the CEP. This effect of the Polyakov 
loop is driven by the fact that the one- and two-quark Boltzmann factors are controlled by 
a factor proportional to $: at small temperature $ ~ results in a suppression of these 
contributions (see Eq. E]) leading to a partial restoration of the color symmetry. Indeed, 
the fact that only the 3— quarks Boltzmann factors e 3/3Ep contribute to the thermodynam- 
ical potential at low temperature, may be interpreted as the production of a thermal bath 
containing only colorless 3-quarks contributions. When the temperature increases, $ goes 
quickly to 1 (this is faster in case I due to the higher momentum quarks present in the 
system) resulting in a (partial) restoration of the chiral symmetry occurring in a shorter 
temperature range. The crossover taking place in a smaller T range can be interpreted as a 
crossover transition closest to a second order one, a feature that is more clear in case I than 
in case II. This "faster" crossover may explain the elongation of the critical region in case I, 
compared to case II, giving raise to a greater correlation length even far from the CEP. 

Now, we will investigate the behavior of Xq an d C in the vicinity of the CEP and their 
critical exponents, for both cases. The calculated critical exponents at the CEP and the 
TCP, together with the universality/mean- field predictions, are presented in Table HVl and 
will be discussed in the sequel. 

In the left panel of Fig. |6] the baryon number susceptibility is plotted as a function of 
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FIG. 6. Left panel: Baryon number susceptibility for case I as function of \i for different temper- 
atures around the CEP: T CEP = 172.48 MeV and T = T CEP ± 10 MeV. Right panel: Baryon 
number susceptibility for case II as function of fi for different temperatures around the CEP: 
T CEP = 169.11 MeV and T = T CEP ± 10 MeV. 

for three different temperatures around the CEP in the context of case I. The behavior 
is very similar in both cases. For T < T CEP , the phase transition is first order and Xg nas 
a discontinuity; for T = T CEP , the slope of the baryon number density tends to infinity 
at ji = fi CEP and Xg diverges; for T > T CEP , the discontinuity of Xg disappears at the 
transition line. A similar behavior of Xg is found for case II, as we can see from the right 
panel of Fig. [61 

The behavior of the specific heat for both function of temperature for three 

different chemical potentials around the CEP, is presented in Fig. [7J The behavior found 
for C around the CEP is very similar to the behavior of Xq f° r both Cases as we can see 
from Fig. [7J 

It is interesting to notice that the high momentum quarks introduced in case I (A — > oo), 
and that are not taken into account in case II, have no significant effect on both e and a. 
However, we observe that the peak at the critical points T CEP or fi CEP is sharper in the 
PNJL model in case I (as it can be expected from the analysis of the stretching of the critical 
region done above). 

To better understand the extreme behavior of Xg an d C near the CEP, we will deter- 
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mine the critical exponents (in our case e and a are the critical exponents of \ q and C, 
respectively). These critical exponents will be determined by following two directions, 
temperature-like and magnetic-field-like, in the (T, /i)-plane near the CEP, because, as 



pointed out in 



45], the form of the divergence depends on the route that is chosen to 



approach the CEP. 

Starting with the baryon number susceptibility, for both cases, if the path chosen is 
asymptotically parallel to the first order transition line at the CEP, the divergence of Xq 
scales with an exponent j g . In the mean-field approximation it is expected to find 7 g = 1 for 
this path. For directions not parallel to the tangent line, the divergence scales as e = 2/3. 
These values are responsible for the elongation of the critical region, Xq being enhanced in 
the direction parallel to the first order transition line (see Fig. 

To study the critical exponents for the baryon number susceptibility (Eq. [T8|) we will 
start with a path parallel to the /i-axis in the (T, /i)-plane, from lower fi towards the CEP 
(T CEP , fi CEP ). Using a linear logarithmic fit lnx q = — eln |/x — \x \ + c\, where the term 
ci is independent of /i, we obtain e = 0.66 ± 0.01, which is consistent with the mean-field 
theory prediction, e = 2/3. This value is also similar to the value found in case II as we can 
see from Table HVl (see also Ref. j^). 

We also study the baryon number susceptibility from higher fi towards the critical [i CEP 
for both cases. The logarithmic fit used now is In Xq — — e' In |/i — fi CEP \ + c[. Our result for 
case I shows that e' = 0.69±0.01~e. This means that the size of the region we observe is 
approximately the same independently of the direction we choose for the path parallel to the 
/x-axis. Once again, this value is similar to the critical exponent for case II and is consistent 
with the mean- field theoretical prediction e = 2/3 (see the results presented in Table HV|) . 

On the other hand, in the chiral limit (where the CEP becomes a TCP), it is found that 
the critical exponent for Xq, denoted by 7 9 , has the value 7 g = 0.53 ±0.02 (7, = 0.51 ±0.01), 
for case I (II). Again, these results are in agreement with the mean-field value (7, = 1/2). 

Let us pay attention to the specific heat around the CEP. We can calculate the critical 
exponent using a path parallel to the T-axis in the (T, /x)-plane from lower T towards the 



CEP. We observe that, as already found in Ref. [27| for case II, the slope of data points 
change for \T — T CEP \ around 0.3 MeV. So, for case I (II) we obtain the critical exponents 
a = 0.64 ± 0.01 (0.62 ± 0.02) and a y = 0.58 ± 0.01 (0.52 ± 0.01). As pointed out in (lo|, this 
change of the exponent can be interpreted as a crossover of different universality classes, 
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FIG. 7. Left panel: Specific heat for case I as a function of T for different values of \i around the 
CEP: n CEP = 286.35 MeV and fx = fi CEP ± 10 MeV. Right panel: Specific heat for case II as a 
function of T for different values of \i around the CEP: /x = 321.32 MeV and \x = fi CEP ± 10 
MeV. 

with the CEP being affected by the TCP due to the small physical quark masses. 

The value of a in both cases is consistent with the one suggested by universality arguments 



m [10|: it is expected that Xq an d C should be essentially the same near the TCP and the 
CEP which implies a = e = 2/3. In addition, if we compare ot\ for both cases, we conclude 
that ot\ in case I is closer to a than in case II. It seems that the high momentum quarks 
allowed in case I affect the crossover of different universality classes making a\ getting closer 
to a which is already consistent with a = e = 2/3 

When the critical point is approached from above the following exponents are obtained: 
a' = 0.68 ± 0.01 (0.68 ± 0.01) for case I (II). 

Finally, concerning the behavior of the specific heat around the TCP, we find as shown 
in Table IIVI a = 0.50 ± 0.01 for case I and a = 0.47 ± 0.02 for case II. These values are in 
agreement with the respective mean-field value (a = 1/2). 



VIII. CONCLUSIONS 



We have considered the PNJL model as one of the prototype models of dynamical sym- 
metry breaking of QCD (both chiral and "color" symmetry) and investigated the phase 
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structure at finite T and \i. Evaluating the thermodynamical potential we find the critical 
curves on the (T, /i)-plane. Working out of the chiral limit, a CEP which separates the 
first and the crossover line is found. First and second derivatives of the thermodynamical 
potential are also evaluated. 

The critical phenomena related with the explicit breaking and partial restoration of Z 3 
and chiral symmetry are analyzed. Two different sets of parameters have been used with the 
emphasis on the parameter choice which is compatible with the formation of stable droplets 
at zero temperature. The effects of two regularization procedures at finite temperature, one 
that allows high momentum quark states to be present (I) and the other not (II), have also 
been discussed. 

We have found that the presence of the Polyakov loop provides a substantial enhancement 
of the critical temperature, bringing it to a better agreement with the most recent results 
of lattice calculations. Another interesting effect of the coupling of quarks to the Polyakov 
loop is that the phase transition becomes steeper, showing a sharper peak in the baryon 
number susceptibility and the specific heat. This effect is reinforced when regularization I 
is used. 

The observation of differential observables, like the entropy and the heat capacity and 
its temperature and density behavior, serve as important probes that, together with lattice 
data, are important for the phenomenology of heavy-ion collisions and cosmology. In both 
cases A and B, as well as in both cutoff procedures, the gross structure of the phase diagram 
expected for QCD is obtained. The location of the CEP depends on the model parameters 
and, in the chiral limit, a TCP is found according to universality arguments. 

Two important points of our model calculation concern the choice of the model parameters 
and the regularization procedure at finite temperature as already referred. We conclude that 
the choice of the model parameters has important consequences in order to obtain the correct 
asymptotic low temperature behavior. In the zero temperature limit, the chemical potential 
approaches a finite value that must satisfy to the condition fi c < M vac . Only the set of 
parameters A insures this condition that allows us to obtain both s = and p q = 0. The 
regularization procedure is important for obtaining agreement with the asymptotic behavior 
above T c . 

Observables like x an d C, which are obtained as derivatives of the thermodynamical 
potential with respect to fi and T respectively, allow us to explore the effects of the Polyakov 



24 



loop on the thermodynamical properties. The successful comparison with lattice results 
shows that the model calculation provides a convenient tool to obtain information for systems 
from zero to nonzero chemical potential which is of particular importance for the knowledge 
of the equation of state of hot matter, relevant for the upcoming LHC experiments at CERN, 
and for dense matter, relevant for the CBM one at FAIR. 
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